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<^> We study the sea quark contribution to the BFKL kernel in the framework of Mueller's dipole 

model using the results of our earlier calculation. We first obtain the BFKL equation with 
the running coupling constant. We observe that the "triumvirate" structure of the running 
coupling found previously for non-linear evolution equations is preserved for the BFKL equation. 
In fact, we rederive the equation conjectured by Levin and by Braun, albeit for the unintegrated 
gluon distribution with a slightly unconventional normalization. We obtain the leading-Aj 
contribution to the NLO BFKL kernel in transverse momentum space and use it to calculate 
the leading- iVy contribution to the NLO BFKL pomeron intercept for the unintegrated gluon 
t>- ■ distribution. Our result agrees with the well-known results of Camici and Ciafaloni and of 

Fadin and Lipatov. We show how to translate this intercept to the case of the quark dipole 

t— i \ scattering amplitude and find that it maps onto the expression found by Balitsky. 

\Q . 

O ■ 
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§^ 1 Introduction 

>• . Understanding the NLO corrections to the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [1,2], the 
Balitsky-Kovchegov (BK) [3-7] and the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov- 
Kovner (JIMWLK) [8-15] evolution equations is important to have their numerical predictions 
under control. Improving the precision of the solutions of BK and JIMWLK evolution equations 
has far-reaching consequences for our understanding of the QCD physics at HERA, RHIC, 
LHC colliders and at the proposed eRHIC/EIC collider. It would sharpen the predictions of 
the saturation/Color Glass Condensate (CGC) physics [3-26] for such observables as total and 
diffractive cross sections, particle yields and spectra, and particle correlations. This would move 
the CGC physics from the qualitative phase to the domain of quantitative predictions. 

Recently there has been much progress in calculating the running coupling corrections for the 
JIMWLK [8-15] and the BK [3-7] evolution equations [27-29]. The idea behind the calculation 
of running coupling corrections in [27,29] was to calculate the quark (leading- Nf) contribution 
to the running of the coupling and recover the leading-order QCD beta-function by replacing 
Nf — > —6 71/32 hi the result. ($2 is given by Eq. (1381 below.) 

The leading- Nf contribution to the running of the coupling was calculated in [27, 29] by 
inserting quark loops to all orders in the gluon line emitted per one step of the JIMWLK or 
BK evolution. If one limits the expansion to just one quark loop, one would then obtain the 
leading- Nf contribution to the NLO JIMWLK and BK kernels. Since the linearized versions 
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of the JIMWLK and BK equations give the BFKL equation [1,2], the results of [27,29] contain 
all the ingredients needed to calculate the leading-iVj NLO correction to the BFKL pomeron 
intercept. To obtain the leading-iV/ NLO BFKL equation from the NLO JIMWLK and BK 
equations with running coupling corrections one should simply keep only the linear terms in 
the latter equations. 

It is a well-known problem that the NLO correction to the BFKL pomeron intercept cal- 
culated in [30-34] is large and negative (for other problems of NLO BFKL see [35,36]). The 
largeness of the correction was argued to be due to collinear singularities [37,38]. Since the 
gluon saturation effects in both the BK and JIMWLK evolutions cut off the contributions from 
the dangerous infrared region, there is a hope that NLO corrections to the BK and JIMWLK 
equations would be numerically small. The work presented below is the first step aimed at 
verification of the above hypothesis in the sense that it explicitly demonstrates the relationship 
between the JIMWLK/BK and the BFKL equations beyond leading order. 

In [29] the leading- iVy contribution to the NLO correction to the BFKL pomeron intercept 
was calculated for the scattering amplitude of a dipole on a target. The calculation was per- 
formed in transverse coordinate space leading to an intercept different from the one previously 
obtained by Camici and Ciafaloni in [34] and by Fadin and Lipatov in [30] using conventional 
momentum space perturbation theory. Below we address the origin of this discrepancy by cal- 
culating the leading- Nf contribution to the NLO BFKL intercept for a different observable: we 
will deal with the unintegrated gluon distribution, defined below by Eqs. (fT5|) and dSJ). We 
show that the leading-iVy contribution to the NLO BFKL intercept for the unintegrated gluon 
distribution obtained below in fact agrees with the results of [30,34]! 

In [27] the running coupling corrections to the JIMWLK/BK evolution equations were 
calculated to all orders. This allows us to derive the BFKL equation including the running 
coupling corrections. We will do this below by keeping only the linear term in the running 
coupling JIMWLK/BK equations obtained in [27]. The result, shown in Eq. fT46|) . gives a first- 
ever derivation of the BFKL equation with the running coupling corrections resummed to all 
orders. While Eq. (1461) is written for the unintegrated gluon distribution, a slight change in 
the normalization of that quantity shown in Eq. (|47p leads to Eq. (1481) . Eq. (T4"8"l) is exactly the 
equation which was conjectured by Braun and by Levin in [39, 40] by ad hoc postulating the 
bootstrap condition to remain valid even when the running coupling corrections are included. 
In this work we present a first-principle derivation of the equation conjectured by Braun and by 
Levin, confirming its validity for a quantity slightly different from the conventional unintegrated 
gluon distribution. A partial validity check of Eq. (1481) was previously performed in [41,42]: 
there the leading-iVy contribution to the NLO BFKL intercept was calculated after expanding 
Eq. (T48]) to the next-to-leading order and was found to be exactly the same as obtained in [30,34]. 

The paper is structured as follows. In Sect. [2] we establish the calculational framework with 
a leading order comparison of JIMWLK/BK on the one hand and BFKL on the other hand. 
We start with the Mueller's dipole model [43-46] analogue of the BFKL equation. We present 
the subtleties involved in performing the Fourier-transform in into transverse momentum space 
and derive the resulting equation for the unintegrated gluon distribution function. We show 
that this evolution equation is indeed equivalent to the standard BFKL equation by calculating 
the eigenf unctions and eigenvalues of the kernel. 

In Sect. [3]we use the results of [27] to construct the leading- Nf NLO correction to the BFKL 
kernel for the evolution of the unintegrated gluon distribution in the large- iV c approximation. 
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The evolution equation combining the LO BFKL kernel and leading- iVj NLO correction to the 
kernel is given in Eq. (I3"9"j) . 

In Sect. H]we construct the BFKL equation with the running coupling corrections included 
to all orders. The resulting equation for the unintegrated gluon distribution is shown in (|4"6]) . 
This is the first main result of our paper. We show that by a simple substitution Eq. (146)) can 
be recast into a form given by Eq. (j4*8l) . which was conjectured by Levin [40] and by Braun [39], 
who postulated the validity of bootstrap condition beyond the leading order. We compare our 
results to the results of [39,40] in Sect. [5] for the case of non-forward BFKL exchange. Again 
we show that we can reproduce the equation conjectured in [39, 40] and clarify the physical 
quantity for which that equation was written. 

We proceed in Sect. [6] by calculating the leading- iVy large- N c correction to the BFKL 
pomeron intercept, which is defined by acting with the leading- iV/- NLO BFKL kernel on the 
eigenfunctions of the LO BFKL kernel. Our result is given in Eq. ( 1681) and is in complete 
agreement with the results of Camici and Ciafaloni [34] and Fadin and Lipatov [30]. We have 
thus verified the results of [30,34] using an independent method. This is the second main result 
of our paper. 

In Sect. [7] we demonstrate that the difference between our result and that of Balitsky [29] 
is solely due to the fact that we are calculating evolution of different observables: while in this 
paper we have been dealing with the unintegrated gluon distribution, the intercept of [29] was 
calculated for the dipole cross section. This difference demonstrates how different the NLO 
BFKL intercept can be for different physical observables! 



2 The dipole model BFKL equation in transverse mo- 
mentum space 

2.1 Fourier transform 

The conventional dipole model analogue of the BFKL equation [1,2] reads [43-46] 
dN{x Q ,x l ,Y) a s N c f 2 x 2 0l 

t^tt = -rr^r dx 2— — r [N{x ,X2,Y) + N[x 2 ,x 1 ,Y) - JV(x ,a;i,y)] , (1) 

0i l, 7T J *^20 "^21 

where x^ = x,i — Xj and = \Xij\. First, for simplicity, we assume that the forward scattering 
amplitude of a dipole on a nucleus does not depend on the impact parameter of the dipole and 
is also independent of the dipole's orientation in the transverse plane: 

N(x Q ,x h Y) = N(x 0h Y). (2) 

The information contained in this simplified object is sufficient to compare to the forward, 
angular-averaged BFKL equation in momentum space. 

To transform Eq. ((TJ) to transverse momentum space let us define 

N(x 01 ,Y) = J-0p (1-e^' 1 ) N(k,Y). (3) 
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Indeed this is not the only way to transform Eq. ([T]) into momentum space: see [7] for an 
alternative. However, it appears that Eq. has the most straightforward application for the 
NLO kernel. 

We recall that N has the interpretation of a dipole scattering cross section, and as such, 
must vanish in the local limit: lim N(xq, x\, Y) = 0; zero size dipoles do not interact. 

This can not be achieved by the using just a naive Fourier transformation, say via the second 
term in ([3]) alone. In fact, if one expands the generic N at low densities, one obtains 



N(x , Xl ,Y) 




I 



Ft 



higher orders 
(4) 



where the horizontal lines at Xq and X\ serve to mark the transverse positions of the quark and 
antiquark. (Note that the complete contribution in (j4j) is positive as it is an absolute value of 
the difference of two amplitudes squared. Similarly the two first diagrams are positive.) The 
blobs in the diagrams are related to the two point correlator of the target fields (integrated over 
longitudinal positions). The relative sign between the two types of contributions ensures that 
the local limit vanishes. This property is preserved if we Fourier-transform the 4 contributions 
individually. To this end we identify (thereby restricting ourselves to the forward case) 



k | — ft 

N(k,Y)= i:i 



(5) 



and take Fourier-transforms of each term in (jlj) with respect to the relative momentum k 
according to 



£ 1 £ 

I i 



d 2 k 

{2-Kf 



k £o i fo —k 
jo yo 



e %h-{xi-Xj) jo I jo 



(6) 



The first two diagrams in (JJ|), with their degenerate Fourier factors e lk ' x °° and e lk ' Xl1 contribute 
the +1 term. For the second term we also use that N(k, Y) only depends on k = |fe|, following 
simplifying assumption made at the outset of this section. 
Substituting Eq. (j3J) into Eq. (TXT) yields 



d 2 k 



Jk 



. xai , dN(k,Y) _ a s N c 



x / d 2 x 2 



x 



01 



x 20 x 21 



dY 
[1 



2tt 2 



d 2 k 



N(k,Y) 



-ikx\2 I ik-XQi 



]■ 
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To obtain a momentum space expression we now use ((7|), along with the Fourier-representation 
of the dipole kernel 



■r, 



01 



X 20 X 21 



d 2 qd 2 q' 



A" 



a. 



g q 

q 2 q' 2 
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and insert both into the low density expanded BK equation ([I]) to obtain 

d 2 k , pik . X01 , dN(k,Y) N c f f <Pqd\' <fk 

\ l ~ e ) m7 = / d x * / ,0^2 (0^2 K i>i' 



(2tt) 2 V 1 dY " 2tt 2 J 2 J (2tt) 2 (2tt) 2 9,9 

x ^iq xo _ giq-asi^ ^g-iq'-asi _ g-iq'-xo^ e -i(q-q / )-x 2 



X 



[(1 - e ik - X02 ) + (1 - e - ifc yi2 ) - (1 - e lk - X01 )] N(k, Y) . (9) 



The goal now is to write the right hand side in a form that resembles the structure of the left 
hand side: 

/d 2 k 
- — rr (l — e tk ' Xtn ) (kernels and N dependence of r.h.s.) . 
(27T) 2 

Once this is done, one may read off the momentum space equation for N(k, Y) by equating the 
integrands. 

To do so one first integrates over X2 to eliminate q' and then transforms the remaining 
contributions (by shifting q and reflecting its sign as needed on a term by term basis) until the 
result for the right hand side takes the form 



2tt 2 



J d 2 qd 2 k { (1 - e tm ) K ieal (q, k) N(k, Y) + (l - e lk ^) K vilt (k, q) N{k, Y)} . 



Last, we exchange integration variables k <h-s> q in the real term and equate the integrands of 
the resulting equation: 



dN(k, Y) N c 



dY 2tt 2 



J d 2 q jiT rea i(fc, q) N(q, Y) + K virt (k, q) N(k, F) j . (10) 



Without using any specific knowledge about the form of the kernel K qq we find that the Fourier- 
kernels K q q only enter in combinations 

L qA :=K qq -K qA (11) 

that vanish in the limit q — > q. In terms of L qq the kernels take the form 

K Tea _\(k, q) =L q ^ k _ k + L_fc,q-fc + L kk ^ q + L k ^ qk (12a) 
K virt (k, q) =L q ^ k + L q , q+k . (12b) 

We will exploit this generic structure below when we consider running coupling corrections. 
For now, substituting the leading order form for K qq into fTTTTl yields: 



a s f 1 1 (q-q) 2 , 



This then is used in (llip and ffTUl) . It is worth noting that only the part symmetric under 
exchange of q and q contributes to ffTU]) . the other contributions cancel out. We find 



dN(k, Y) a s N c 



d 2 q 



la 2 k 2 

N(q,Y)-—- ^N(k,Y) 



8Y 2tt 2 J * [k 2 (k-q) 2 w ' q 2 (k-q) 2 



(14) 
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Finally, defining the unintegrated gluon distribution via 



AT Q 

a s (k 2 ) <f>(k,Y) = -J—Lk 2 N(k,Y) (15) 



we write 

d<p(k,Y) a s N, 



d 2 q 



<Kq,Y)- 9/I fc2 , 9 0(fc,y) 



dY 2tt 2 J * L( fc -9) q 2 (k-q 



(16) 



At the leading logarithmic order we have neglected the running of a s and dropped this factor 
on both sides of Eq. ([T6"j) . Eq. (1161) is the momentum space representation of the dipole model 
analogue of the BFKL equation. 

The definition of the unintegrated gluon distribution <p(k, Y) given by Eqs. (I15p and ([3]) 
is illustrated by the figure of Eq. (j3j). Since we are interested in the linear (low color density) 
regime the interaction between the dipole and the target is limited to the exchange of a single 
BFKL (or, as we will consider later, NLO BFKL) ladder. Thus the unintegrated gluon dis- 
tribution <p(k, Y) includes the whole (NLO) BFKL ladder attached to the target. Summation 
over all possible connections of the t-channel gluons to the quark and anti-quark lines in the 
dipole has to be performed, as shown in Eq. (J4j). Due to our assumption of independence of 
the dipole amplitude N on the impact parameter shown in Eq. ([2]) we were able to explicitly 
integrate over the impact parameter in Eq. (TIB]) which resulted in the factor of the transverse 
area of the target S±. 

The coupling constant a s (k 2 ) accounts for the interaction of the t-channel gluons with 
the quark lines of the dipole: indeed this interaction is not a part of the unintegrated gluon 
distribution function and has to be factored out. The scale of the running coupling is naturally 
given by k 2 as the only available scale in the problem. As we will see later, even at the NLO 
BFKL level one does not need to know the constant under the logarithm in this factor of 
the running coupling. If one takes the diagrams in Eq. (Jlj) literally and tries to calculate the 
contribution of the two exchanged gluons with the running coupling corrections included, one 
would get a factor of (a s (k 2 ) 2 /a M ) l/k 4 , with a M the bare coupling constant. Here we follow 
the standard convention for the unintegrated gluon distribution function and require that at 
the lowest (two gluon) level it should be given by 

Mk,Y) = ^' (17) 
7i k 

Unintegrated gluon distribution defined by Eq. (IT?]) actually gives us the number of gluon 
quanta in the phase space region specified by its arguments. To adhere to this lowest order def- 
inition of the unintegrated gluon distribution we absorb the factor of (a s (/c 2 )/a M ) 1/k 2 coming 
from the two gluons in the figure of Eq. (jlj) into <fr(k, Y), leaving the factor of a s (k 2 )/k 2 out in 
front, as we see in Eq. ([TBI) . 



2.2 Leading order intercept 

Eq. (fT6l) is equivalent to the BFKL equation [1,2], averaged over azimuthal angles. To see that 
this is indeed the case, and to prepare for the pomeron intercept calculations to be carried 
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out below, let us show that powers of momentum are indeed the eigenfunctions of the kernel 
of Eq. (1161) and find the corresponding eigenvalues. Acting on (p(k, Y) ~ k 2X with the kernel of 
Eq. (US} yields 



a x N r 



7T 



J d 2 qK LO (k,q)q 



x _ a s N c 
2n 2 



d 2 q 



,2 A 



k 2 



„2A 



(k — q) 2 q 2 (k — q) 



Rewriting the measure of the q-integral as [43] 

oo 

d 2 q = 2ixdqqdll J db b J (b k) J Q (b q) J (b I) 
'o 

with k = \k\, q = \q\, I = \k — q\ we obtain 

oo oo oo 

J d 2 qK LO (k,q)q 2X = 2 n J dbbJ (bk) J dqqJ {bq) J dll~ l J {bl) 



(19) 



-2A+2 



2q 2 



(20) 



where we introduced an infrared cutoff A to regulate q- and /-integrals. Using Eqs. (IAip and 
( 1A2I) in Appendix [A] we can perform the q- and /-integrals recasting Eq. (1201) into 



J d 2 qK LO (k,q) q 2X 

oo 

= 2 j dbbJ (bk) ^(l)-ln^ 



22A+1 ^-2 A-2 



r(i + A) k 2X+2 
r(-A) 2 



^(l)-ln^ 



(21) 



After performing the integration over b (again with the help of the formulae from Appendix [A]) 
Eq. (T2~T]) reduces to 



a s N, 



71 



c J d 2 qK LO (k iq )q 



a s N c 



7T 



X(-A)fc 



2A 



where 



X (-A) = 2V(1)-V(-A)-V(1 + A). 



(22) 



(23) 



This accomplishes the proof that, at least in the azimuthally symmetric case, Eq. (Tl6|) is 
equivalent to the BFKL equation. 



3 Momentum space NLO BFKL equation with 
the leading- Nf NLO correction in the kernel 

To calculate the leading- Nf contribution to the NLO BFKL kernel in the language of Mueller's 
dipole model [43-45] we have to insert a single quark bubble correction in the gluon line emitted 
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in one step of the dipole evolution. (For a brief review of Mueller's dipole model see [26].) The 
relevant diagrams are shown in Fig. [Hand have been calculated recently by the authors in [27] 
and by Balitsky in [29] 




ABC 



Figure 1: Diagrams contributing a factor Nf at order a 2 ln(l/:r). The line from bottom left 
to top right indicates the target field crossed by the fast moving constituents of the projectile. 
Their transverse positions are fixed. The ovals indicate insertions of N(xq, Z\,Y), N(xq, z, Y) 
and N(x , x±, Y) respectively. 



The calculation in [27] followed the rules of the light cone perturbation theory (LCPT) 
[47,48]. The diagrams in Fig. [T]should also be understood as LCPT diagrams. The propagators 
of gluon lines in the graphs A and C in Fig. [1] also include the instantaneous parts. 

Here we are interested in the linear evolution: hence only one of the produced dipoles in 
Figs. [T] A and B continues to have evolution corrections and interacts with the target. In Fig. [T] 
we show that with an oval denoting the interacting dipole. Indeed either one of the dipoles can 
interact and later on we will add the contribution with the other dipole interacting with the 
target. 

We start with the contribution to the NLO JIMWLK kernel coming from the diagram A in 
Fig. [1], which is given by Eq. (17) in [27]: 



K^ hO (x ,x 1 ;z 1 ,z 2 ) 



+ 



da 



d 2 k d 2 k' d 2 q d 2 q' 
(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 " 

2a) 2 q k k' q' + q q' k k' 



-iq-(z—x )+iq'-(z—xi)—i(k—k')-zi2 



k' k 



q 2 q' 2 



k 2 

2a(l - 



q 2 a(l 



a) 



a 



2a) 



k' 2 + q' 2 a{l - 
k ■ q 



a 



k 2 + q 2 a(l 



a 



k 



12 



q' 2 a{l 



a) 



+ 



k'q' 



t'2 



1 We emphasize that the derivation there was performed in the JIMWLK/BK context that includes nonlinear 
contributions. Here we only retain the linear terms relevant for the BFKL limit. Using light-cone perturbation 
theory one might have calculated these contributions directly in the dipole model. 
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4 a 2 (1 -a) 2 

+ 7 TT — 

k 2 + q 2 a(l — a) k' 2 + q' 2 a(l — a) 



(24) 



At the moment we exclude all the factors of the coupling constant and the color factor, similar 
to what was done in [27]. The notations used in Eq. (j2"4"|) are explained in Fig.[I]A the quark and 
the anti-quark have transverse momenta ki and k 2 and transverse coordinates Z\ and z 2 . The 
gluon line carries momentum q in the amplitude. All the transverse momenta in the complex 
conjugate amplitude are labeled in the same way as momenta in the amplitude only with a 
prime. The longitudinal fraction of the gluons' momentum carried by the quark is labeled 
a = k\ + /{ki + + k 2+ ). We also note that z = a z\ + (1 — a) z 2 , Z\ 2 = Z\ — z 2 , q = k± + k 2 and 
k = fex(l — a) — k 2 a. 

Using q = ki + k 2 and k = k\(l — a) — k 2 a (and the same for primed ones) we rewrite 
Eq. ([MD as 

l 

P 1" d 2 k, rl 2 k n rl 2 h' rl 2 h' 

NLO / . \ A AT J I 1 2 1 2 



x e 



V? uu (x ,x 1 ;z 1 ,z 2 ) = ANf J da 

o 

-iki-(z 1 -x )—ik2-{z2-x )+ik' 1 -(z 1 —x 1 )+ik' 2 ■(z 2 ~x 1 



(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 
1 



(k, + k 2 ) 2 (k[ + k' 2 ) 2 

k x ■ fei k 2 -k' 2 - fci • fc^ fci • k 2 + [fci ■ k 2 + fc 2 (1 - a) + fc| a] [fci • k' 2 + fcf (1 - a) + fc^ 2 a] 



[k 2 (1 - a) + k 2 a] [kf (1 - a) + & 2 2 a] 



(25) 



As indicated in Fig. [TA, in the linear regime only one of the dipoles, say the dipole formed 
by Xq and Z\ continues the evolution. Hence there is no dependence on the coordinate z 2 in the 
resulting iV(|:Eo — zjj, Y). Therefore we can (and should) integrate z 2 out. Adding the result 
to the other Nf corrections leads us to consider the whole contribution as part of the running 
coupling corrections. 

Let us note that this procedure is equivalent to the "UV subtraction" performed by Balitsky 
in [29] for the nonlinear evolution equation and briefly explain the context. For the nonlinear 
equations, extraction of the UV-divergent part of the diagram A in Fig. HJ which one must 
absorb in the running coupling constant, can be done in many ways, depending on which linear 
combination of Z\ and z 2 is kept fixed. This choice is referred to in [27], maybe somewhat 
unfortunately, as the choice of "subtraction point." Since the notions of "UV subtraction" 
and "subtraction point" play an important role in the renormalization of QCD, we should 
attempt to forestall unnecessary confusion by pointing out that the ideas involved here are only 
peripherally connected to renormalization. The subtractions referred to here correspond to 
different schemes of separating the running coupling contributions (which must carry the UV 
divergence of QCD diagrams so that the UV scale \x will enter the logarithms that induce the 
running of the coupling) from new physics channels (such as, in this case, the presence of a well 
separated gg-pair in the final state), from which the UV divergent running coupling contribution 
has been subtracted. The only new channel we have considered explicitly, the qq channel, then 
turns out to be completely finite and \i independent at present accuracy. This separation 
scheme dependence (and the coordinate space "subtraction points" associated with it) is a 
freedom that arises independently in addition to the renormalization scheme dependence and 
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the momentum space renormalization scale /i (with associated subtraction points) appearing in 
the renormalization procedure of QCD. In the case of non-linear evolution the arbitrariness of 
the separation scheme results in the differences between [27] and [29]. We emphasize that the 
total result, in which running coupling contributions and new channels are added together, is 
not affected by this separation. In contrast to renormalization scheme dependence, which can 
not be eliminated at any finite order in perturbation theory, this statement about separation 
scheme independence holds order by order in perturbation theory. In the separation of [27] both 
the running coupling corrections and the new qq channel contribute leading-iVy terms to the 
linear part of the BK equation and the sum of these contributions is obtained by the the above 
procedure. In the separation of [29], all leading order Nf- contributions to the linear part of the 
BK equation are included into the running coupling, while the correspondingly constructed new 
qq channel only contributes a nonlinear term to the BK equation. This establishes the direct 
correspondence of our present procedure (which is simply aimed at collecting all leading- Nf 
contributions to the linear part of the BK equation) to the subtraction used in [29] . 

Let us stress again that for the total leading- Nf contribution (and that is what we are 
calculating here) subtraction points are not an issue, they merely affect the separation into 
what we interpret as a running coupling correction and new production channels. The only 
remarkable point here is that the leading- Nf contributions may all be interpreted as running 
coupling corrections. 

To return to our present task: we are interested in the linear evolution and have to integrate 
the kernel in Eq. (125]) either over z 2 or Z\ depending on which dipole interacts with the target. 

To integrate over z 2 , which would put k' 2 = k 2 , we define q = k\ + k 2 and q' = k[ + k' 2 = 
k[ + k 2 , obtaining 

l 

J d z 2 k, l (x ,x 1 ,z 1 ,z 2 ) - 4 Nf J da J e 



q 2 q' 2 [(q - k) 2 (1 - a) + k 2 a] [(q 1 - k) 2 (1 - a) + k 2 a] 
x <j fc 2 (q-k) ■ {q' - k) - k- (q- k) k ■ (q' - k) + [k ■ (q - k) + (q - kf (1 - a) + k 2 a] 

x [k ■ {q -k) + (q 1 - kf (l-a) + k 2 a 



(26) 

where we dropped the subscript 2 in k 2 , i.e., k 2 — > k. Now we have to integrate over k. To do 
so we first write 

i 

( d 2 Zo JC NLO (xn xv z^ Zo) = 4N f fda I - 9 - q ' - - e -iQ<^-^o)+ig-( Z i- X i) 
7 2 1 WoiXxiZuZz) 4i V/ j aa j (2?r)2 (27r)2 (2?r)rf 



if k 2 {q-k)- (q f - k) 

X ~qW 2 | [{q ~ k) 2 (!-«) + k 2 a] [{q' - k) 2 (1 - a) + k 2 a] 
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k-(q-k) k {q' - k) 

(q-k) 2 (l-a) + k 2 a + {q' - k) 2 (1 - a) + k 2 a 



(27) 



where, in anticipation of dimensional regularization we have inserted the dimension d of the 
fc-integral explicitly. The details of fc-integration in Eq. (p7|) are given in Appendix [Bl The 
result yields 

1 f Co - o'l 2 P" 5 / 3 „2 "5/3 /2 -5/3] 

x ^ <«-«') 'i-iHe — " 2ln ^ — 9 ' 2ini ^ ■ (28) 

q q { ^MS ^MS ^MS J 

Note that due to the symmetry property 

K* hO (x , aj i; Zi, 2 2 ) = /Cf LO (* , xi; «i) (29) 

the contribution to the NLO BFKL kernel of Fig. [TJA. with the dipole formed by z 2 and X\ 
interacting is equal to the right hand side of Eq. (128j) with Z\ — > 22- 

To construct the full leading-iV/ NLO BFKL kernel one has to add to Eq. ( 1281) the contri- 
bution to the NLO BFKL kernel /C^ LO of the diagram B from Fig. [TJ Using Eq. (28) in [27] 
along with its mirror reflection with respect to the line denoting the interaction with the target 
we obtain 



/ 



d 2 z 2 K^ o (x ,x 1 ;z 1 ,z 2 )+]C^ o (x ,x 1 ;z 1 ) = ^ / e -M— 

(9-9) ln 2 9 In 2 q In - 2 

9 2< ?' 2 Rm 



MS HVIS ^MS 

2 -5/3 fl /2 -5/3 

2 q ■ q' ln +2q-q' \n 



(30) 



The above expression can be simplified to give 

72^ j2„i 



y 3 7T J (27T) 2 (27T) 2 



9 In h q In h 2 q • q ln 



^2 g/2 j ^2 g/2 (q — q') 2 n 2 



(31) 

where we have also defined /1 2 = /i^e 5 / 3 . 

We do not need to calculate explicitly the diagram in Fig. [1] C due to the real-virtual 
cancellations which lead to the identity 

J d 2 z 2 IC* LO (x , x x - Zl , z 2 ) + K^°(x , xi; z x ) + K* LO (x , Xl ; z x ) = 0. (32) 
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Here 



J d 2 z 1 JC^°(x ,x 1 ;z 1 ] 



is the contribution of Fig. [T] C represented as the integral over the transverse coordinate of the 
(anti-)quark. 

The full leading-Afy NLO dipole kernel is defined via 



,x 1 ;z 1 ) = y Yl (-l) m+n { / d 2 z 2 K^°{x m ,x n] z 1 ,z 2 )+K^°{x m ,x n - 1 z l )\ 

m,n=0 ^ 



(33) 



where we have included the color factor N c /2 in the large- N c approximation and have summed 
over all possible emissions of the gluon in Fig. [1] off the quark and the anti-quark lines in the 
original dipole. 

Including the contributions of either one of the dipoles in Figs. [T] A and B interacting with 
the target we derive that the additive NLO contribution to the right hand side of Eq. ([T]) is 
given by 

al J d 2 x 2 Kf LO (x , x x] x 2 ) [N(x Q , x 2 , Y) + N(x 2 , x u Y) - N(x Q , x u Y)] , (34) 

where a M is the bare coupling constant. In the discussion of LO BFKL the running of the 
coupling was negligible: hence we did not distinguish the bare coupling from the physical 
coupling and labeled them a s . At the NLO level running coupling corrections become important 
and we will start distinguishing between the two. 

Combining Eqs. (13~T1) and (13"3"|) yields an additive NLO modification of the LO dipole kernel 
from Eq. (jSJ) of the form 

K q,q' <*/i 2 —f 1 • 35 

Since our manipulations leading to the evolution equation (flOl) did not make use of any special 
properties of the kernel, we only need to form linear combinations according to (j!2p and ffTTl) to 
find the evolution equation for the unintegrated gluon distribution function defined in Eq. (FTST) 



2 d(j)(k,Y) a^N c . 2 

a *( k ) ^7 = 9 / d 1 



dY it 2 
d 2 q 



(f)(q, Y) -a s {k) -— <j>{k, Y) 



(k - q ) 2 ' q 2 (k- q] 



cjN c N f 

12 7T 3 



k — q) 2 \ q 2 ji 2 J ' q 2 (k — q) 2 \ k 2 \i 2 



(36) 



Here we could not neglect the running of the coupling and kept the factors of the physical 
coupling a s introduced in the definition of the unintegrated gluon distribution in Eq. (fTol) . 
However, as we are interested in the NLO correction only, we should divide both sides of 



12 



Eq. (1361) by a s (k 2 ) and expand the kernel on the right hand side of the resulting equation up 
to the second order in using 



a^Nf k' 
1 - -g— ^ In - 

7T 0^ 



(37) 



Here 



02 



11 N c - 2 N f 
12tt 



and we replaced — » —Nf/6n since we are interested in leading- iVj contribution only. 
With the expansion of Eq. 037|) we finally obtain 



d<j>(k,Y) _ a^N c 



dY 



x / d 2 q 



(k-q) 



2tt 5 



In 



d 2 q 



I? 



(k-q) 2 

<f>(q,Y) 



cf>(q,Y) 



q 2 (k - q) 



<f>(k,Y) 



+ 



k 1 



q 2 (k - q) 



In 



q 2 (fc ~ g) 
A; 2 /i 2 



12 7T 3 



(3? 



(39) 



Eq. fl3T?|) is the BFKL equation with the leading- Nf NLO correction calculated in the large- N c 
approximation. 



4 Resummation of bubble diagrams: triumvirate 
couplings for forward BFKL 

The leading-TV/ NLO BFKL kernel in Eq. fl39l) is remarkably simple, appearing to be more 
compact than the similar kernel obtained using conventional perturbation theory in [30,34]. To 
understand the origin of this simplicity let us point out that, using the techniques developed 
in [27] one can resum the a M Nf corrections to all orders, which correspond diagrammatically 
to inserting an infinite chain of quark loops on the gluon line emitted in one step of small- x 
evolution in the s-channel approach. Replacing Nf — > —6 tt (3 2 and absorbing all the corrections 
into the running coupling constant we obtain the JIMWLK kernel with resummed running 
coupling corrections □ 



a M /Crc(aj ,aJi;«) = 4 



d 2 q d 2 q' 
(2tt) 2 (2tt) 



-iq-(z-x )+iq' -(z—Xj] 



q ■ q 

q 2 q' 2 



a . 



(q 2 e~ 



5/3 ) a s 



(q 



12 -5/3 



) 



a s (g 2 e" 5 /3) 



(40) 



where 



In 



Q 2 e 



-5/3 



/4s 2< ?<7' 



2 -5/3 12 -5/3 

2 1 q e a l q e 

q m ^ ^ q m ' 



In 



(9 



q') 2 e" 5 / 3 



(41) 



2 We refer to this kernel as the JIMWLK kernel since it can be used to construct the JIMWLK evolution 
equation with the coordinate space subtraction point of [29]. Below we will use its Fourier transform to 
momentum space as the kernel of the BFKL equation with running coupling corrections. 
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This expression leads to a remarkably simple form for the building blocks of the forward kernels 
for the angular averaged case: We find 

NLO _ 1 f a s {q 2 e- 5 / 3 ) a s {q' 2 e^ 3 ) a s (q 2 e-^)a s (q' 2 e^ 3 ) (g - q') 2 
2\ q 2 q' 2 a s ((q-q') 2 e" 5 / 3 ) q 2 q /2 

nlo 1 f «*(<? 2 e" 5 / 3 ) a s (q' 2 e" 5 / 3 ) , a s (q 2 e" 5 / 3 K(q' 2 e" 5 / 3 ) (q - q') 2 



L ™'^ L ™' 2\ q 2 q' 2 a s ((q-q') 2 e~ 5 / 3 ) q 2 q' 2 



(42a) 



(42b) 



Since the kernel in the real part is an explicitly symmetric combination of L's, one can imme- 
diately read off (dropping the NLO superscript for brevity) 

T + T T T , a s (k 2 e~^) a s ((g - k) 2 e^) q 2 

For the virtual contribution one notes that the antisymmetric parts of both terms integrate to 
zero (by a shift in q). We reverse the sign of q in the remainder of the second term and are 
left with 

U 2 a(h + L \ - f d 2 a ^(q 2 e-^)a s ((q~kre-^) k 2 
Jdq^ q _ k + L q>q+k )- Jdq as(k 2 e - 5/3) q*(q-k) 2 ' {U) 

Collecting all the contributions one obtains an equation for N(k, Y): 
dN{k,Y) N c f f a s {ee-^)a s {{q-kfe^) q 2 ~ 



TV^yU. (45) 



<9F 2tt 2 7 * | a s (q 2 e- 5 / 3 ) fc 2 (q-fc) 2 

a s (q 2 e" 5 / 3 )a s ((q - fc) 2 e" 5 / 3 ) fc 2 



a s (fc 2 e -5 / 3 ) q 2 (q — k)' 



Using Eq. (1151) in Eq. ([45]) yields the following evolution equation for the unintegrated gluon 
distribution including running coupling corrections to all orders 



d<f>(k,Y) N c 



dY 2tt 



k 2 a s {q 2 e-^) a s ((fc - q) 2 e~ 5 / 3 ) 
~ q 2 (fc - q) 2 a s (fc 2 e~ 5 / 3 ) 



0(fc,F)L (46) 



Expanding Eq. (1461) in powers of to order a 2 would yield Eq. ( I39i) . Hence, from the stand- 
point of linear evolution, the simplicity of the leading- iV/ NLO BFKL kernel obtained here is 
due to the fact that this kernel consists only of running coupling corrections. 

Eq. (|4"6"1) should be compared with the evolution equation for the unintegrated gluon distri- 
bution derived in [40] (c.f. Eqs. (22) and (23) there). At first sight Eq. (146]) appears to disagree 
with these expressions. However, defining a new function 

= r£-L ( 47 ) 

a s (k z e 
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we recast Eq. ( 146|) into 

d4>(k,Y) N c f 2 a s (q 2 e^ 3 ) a s {{k-qfe^) 
dY ~ 2vr 2 J q a s (fc 2 e- 5 /3) 

^v^^'-^i^w^} (48) 

Eq. ( HHl) is in agreement with Eqs. (22) and (23) of [40]. Eq. (T4T1) links the evolving quantity 
to the unintegrated gluon distribution. 

Eq. (j4"8l) was originally obtained in [40] by calculating running coupling corrections to the 
virtual contributions in the non-forward case. This result was then generalized to include 
running coupling corrections also to the real contribution by postulating the validity of the 
bootstrap equation beyond leading order [39,40] . To understand if this postulate can be justified 
from our calculations, we need to explore also the non-forward case. 



5 Generalization to non- forward BFKL: triumvirate 
couplings and bootstrap ideas 

Now that the general techniques are already familiar we relax our restrictions to forward BFKL 
and angular averaged results. This will lead us to abandon our initial assumption of impact 
parameter independence. Practically this amounts to use a Fourier-representation of the lowest 
order expansions of N(xq, a?i, Y) shown in Eq. (jlj) using two independent momenta. We define 



which are related to each other by the double Fourier-transform 

r d 2 k d 2 k' ik , Xi+ik >. x r _ f d 2 k d 2 l ik . ( )+il .^± , . 

After the second equality sign we have changed integration variables to conjugates of the rel- 
ative and absolute positions in a form most suitable to our calculations below. With this 
definition Eq. is replaced by 

N( Xi , Xj , Y) J- J ( g 2 (e fl - + - - e -*-«+*=^) G k+ i_ k _i_(Y) 

+ higher orders (51) 

in a term by term correspondence with Eq. (J3J). Despite this change compared to the forward 
case, one can still proceed to extract the (non-forward) BFKL equation using a technique that 
completely parallels the steps used in Sec. [2j The only price to pay is additional algebraic effort. 



X; 



i 1 1 



X; 



k 




G 



k,k' 



i 1 1 




(49) 
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In this spirit, we insert our expressions for N(xi, Xj, Y) and a generic kernel of the form (jSj) into 
the linearized BK equation ([T]) as before. We will later use the resummed expressions shown 
in (1421) to obtain explicit results. 

In our calculation, we first perform the ^-integral to eliminate the q' dependence and again 
use shifts and sign reflections in q to cast the result in the form 



d f d 2 k d 2 l f ^i. X0 e ; tol _ eifc . aJ01 + il .£!a+2i _ e _ ife . ffi01+a .£jfl+!Ji » G 



dY J (2tt) 2 (2tt) 2 

d 2 k d 2 l N, 



d 2 q | {e llx ° + e* 1 * 1 - e^i+u-^ 1 - e -im+u-^^ j K Ieal (q, k, I) 
_ ^U-wo + e il Xl _ _ e _ ifc . aol -Hi.-D+-t^ # virt ( fe) g> /) j G fc+ i jfc _ i (Y). 

(52) 

The kernels can again be expressed in terms of L qq i as before in the forward cas^l We find 

K Tca \(k, q, I) = L k _^_ q+k - L fe _| )fe+ » - L 9 _ fe _ fc+ j + L_ k _i_ q _ k - L_ k _^_ k+ i_ + L_ q+k k+ i_, 

(53a) 

JCirt(fc, 9, = £ q> _ fc+ ! +q + ^q,fc+I+q- (53b) 

In these expressions we have already anticipated the swap of fc and qr in the real term to 
facilitate the last step, in which we read off the equation for G k+ i k _i(Y) by equating the 
integrands: 



j d 2 q {K T ^(k,q,l)G q+iq _ L2 (Y) - K vilt (k, q,l)G k+hk ^(Y)} (54) 



dG k+! _ k „i_(Y) _ n c 
dY 

Last, we insert the explicit expressions of Eq. (]42p into fl53|) and repeat the arguments of the 
forward case to again cancel the antisymmetric parts. Anticipating a comparison with Levin 
and Braun [39,40] we introduce the notation 

"(* 2 > : =MW) (55) 

and change notation by shifting both k and q by —1/2. Note that [39,40] do not specify the 
renormalization scheme and hence do not keep track of the factors e _5//3 we take care to include 
here@ 

With these conventions, the final equation reads 



dG k , k _ t (Y) _N C 
dY 2tt 2 



J d 2 q K rcal (fc, q, I) G qiq -i(Y) - [0 ((k - if) + (3 (k 2 )] G fe)fe _j(y) (56) 



3 Note that in ([52]) we may reverse the sign of q in K rea \(q, k, I) on a term by term basis. The form chosen 
here reduces the number of steps needed to arrive at our final result (|56|) . 

4 The factor simply enters the definition of Aqcd and emerges directly from our calculations in dimensional 
regularization. This is entirely independent from the issue of separation scheme dependence as discussed in 
Section [3J 
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where the real kernel takes the explicit form 



„ = . "((«-')') 2 ^ + up up (57) 

r,{(k- If) , (((fc - qf) q (fc 2 ) ,((*- g) 2 ) ,,(fc 2 ) If ((fc-() 2 ) 
and /? (fc 2 ) denotes the gluon trajectory defined as 

g(*l4M; rirw ^ - < 58 > 

We note that all individual contributions are composed of triumvirate structures. In the forward 
limit, this reduces to what we found earlier in (|45j) with iV(fc, Y) « G^^Y). 

We now wish to compare to [39,40], without loosing the connection with the bootstrap 
condition used in their argument. This prevents us from taking the forward limit (which we 
have already shown at the end of Sec. @J. Instead, we identify the quantity used in [39,40] by 
amputating the external, coupling-resummed gluon legs from the dipole scattering amplitude 

Gk t k'(Y): we separate out a factor - ""^^ - according to 



k !g 

Gk,k>(Y) = io 




£> k' 

IS = 1 fc2 ^ 1 p ^ ^MY) ■ (59) 



Generically, such modifications leave the gluon Regge trajectories unaffected and rescale only 
the real kernel. Here, the real kernel is modified to 

K real (k, q, I) =K reSL i(q, fc, I) 



7]((k-l) 2 ) 7](k 2 ) 7] (l< 



rj ((q - I) 2 ) V (((fc - q) 2 ) rj (q 2 ) rj ((fc - q) 2 ) rj (q 2 ) rj ((fc - I) 2 ) 
This turns Eq. (156]) into 
dip k ,k-i(Y) N c 



(60) 



dY 2tt 



J d 2 q K rcal (fc, q, I) Vq ^i(Y) - [(3 ((fc - I) 2 ) + (3 (fc 2 )] (fk,k-i(Y) . (61) 



The structures in this equation correspond directly to (16) of [40] or (4) of [39], which, as 
already mentioned, were based on calculating the running coupling corrections to the virtual 
contributions and postulating a bootstrap condition to extend the result to include corrections 
for the real contributions. 

Our calculation now provides a derivation of this statement in the s-channel light cone 
perturbation theory formalism, previously used to derive the JIMWLK and BK equations, and 
gives a firm interpretation of the objects involved. 
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6 Leading- Nj contribution to the NLO BFKL pomeron 
intercept 

To find the leading- iV/ contribution to the NLO BFKL intercept (in the large- N c limit) we have 
to act on <J)(k, Y) ~ k 2X with the NLO part of the kernel in Eq. ( |39|) . Employing Eq. ( fl9l) again 
and using the same notation as in Sect. 12.21 we write 

J d 2 qK^ LO (k,q)q 2X = 

J2 \ 2 , f (k-q) 2 \ k 2 
aq [(k-qy V V 2 ) q ~q 2 (k-q) 2 



In 



q 2 (fe ~ q) 

k 2 fi 2 



„2A 



2tt J dbbJ (bk) J dqqj (bq) J dll' 1 J (bl) X 
A A 



(4) f x -% ^(^- 2 \ 

fi z J q z \k z 



Integrating over q and I with the help of formulae in Appendix |A] yields 

J d 2 qKf LO (k,q) q 2X 



2tt / dbbJ (bk) (^(l)-ln 



bA 



>2A+2 



-2X-2 r(i + A) 
r(-A) 



21 



^(1) - In 



b fi 



(62) 



-2^+ 2 f^(l)-ln^) (W)-m^ 



Finally, using formulae in Appendix |A] to integrate Eq. (IB"3"j) over 6 we obtain 



Including the prefactor of which we have been omitting above, combining (|64|) with the 

leading order contribution to the intercept (1221) and remembering that /i 2 = „ e 5//3 yields 



(63) 



X 2 (-A) + x(-A) In ^ + ~ V'(-A) - ~ ^'(1 + A) \ . (64) 



A/5 



71 



2 A 



^ ' X(-A) 



7T 



1 - a„ fi 2 In 



12vr 



10 



- X 2 (-A) - V'(-A) + ^'(1 + A) + — x(-A) 



(65) 



where we have replaced Nj — > — 67r/3 2 in front of the logarithm to underline the fact that this 
separation of terms has often been interpreted a separation of running coupling and conformal 
contributions. 

From our earlier considerations in Sec. [3l we now know that such a separation is artificial. 
In fact the discussion there and in Sec. H] as well as Sec. [5] allow us to take the extreme position 
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and assign all contributions listed here to running coupling effects. To facilitate comparison 
with [30,34], we still introduce the notation A(A) for the last term in Eq. (I65p by 



7T 



X(-A) 



OluNc 



From Eq. (165]) we have 



l-OLufc In 



12 7T 3 

k 2 ■ 



MS 



d 2 qK^ LO (k, q) 



12vr 



A(A) 



A(A) 



10 



k) 



2 A 



(66) 



- X 2 (-A) - V'(-A) + ^'(1 + A) + — X (-A). 



(67) 



To compare our result with the result of the NLO calculation of Fadin and Lipatov [30] and 
of Camici and Ciafaloni [34] we note that in [30, 34] instead of the power A the authors used 
the power 7 — 1, such that A = 7 — 1. We thus first rewrite Eq. fl67|) in terms of 7 



A( T ) 



-X 2 {l) - ^'(1 " 7) + ^(7) + -7T Xil) 



(68) 



The intercept of Eq. (1681) exactly agrees with the result of Fadin and Lipatov (see Eqs. (14) 
and (12) in [30]) and with the result of Camici and Ciafaloni (see Eq. (4.6) in [34])! We have 
thus provided an independent cross-check of those earlier results [30,34]. 

It is important to note that the agreement between the intercepts of [30,34] and the one 
in Eq. ( 1681) depends crucially on the fact that both intercepts were calculated for the same 



observable — the unintegrated gluon distribution. We also stress that the factor of a s (k ) in 
the definition of the unintegrated gluon distribution in Eq. (fi"5l) plays an important role in 
obtaining the correct intercept. 



7 Comparison with the NLO intercept for the evolution 
of dipole amplitude calculated in [29] 

Another leading-iVy NLO BFKL intercept we should compare with is due to a transverse 
coordinate space calculation performed by Balitsky [29]. The result of [29], expressed in terms 
of the same power 7 as defined in [30,34] reads (see Eq. (44) in [29]) 

4 10 

A B (7) = X\l) + ^(1 - 7) - ^(7) - x{l) + ir x{l)- (69) 

1 — 7 3 

We will demonstrate below that this expression can be translated into the corresponding 
(uniquely defined) contribution in the results of [30,34] and our earlier result Eq. ( 1681) . by taking 
into account that (i) Ab(7) was obtained using both kernels and eigenfunctions in transverse 
coordinate space and that (ii) Ab(7) is the intercept for a different observable — the forward 
scattering amplitude for a dipole on the target, while A (7) was calculated for the unintegrated 
gluon distribution in momentum space. 
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We will, therefore, show that after a Fourier transform relating <p(k,Y) to N(xqi,Y), our 
intercept from Eq. (1681) is translates into the intercept Ab(7) from Eq. (1691) obtained in [29]. 
Combining Eq. with Eq. we write 



27T 

N~c 



k 2 



(70) 



Therefore, if we put a s (k ) (f)(k, Y) = k then we would obtain 



N(x 01 ,Y) 



(27T) 



,2A-1 „-2A 



r(A) 



n c - ^ 01 r(i-A)' 

Using Eq. (|36|) and the steps which led to it we write 

al J d 2 x 2 Kf LO (x ,x 1 ;x 2 ) [N(x Q , x 2 , Y) + N(x 2 , x u Y) - N(x , x h Y)] 



(71) 



2 n fd 2 k 



Nr. 



k 2 



alN c N 



I 



12 7T 3 



d 2 g 



i^ LO (fc,q) + 



In 



k 



21 



(k — q) 2 q 2 



a s (q 2 )^q,Y). 



(72) 



2 A 



Therefore, for N(x i,Y) = x { 



«^r°^oi 2A 



-2A+1 



2vr 



r(l-A) f d 2 k 



k 2 



<xjN e N f 

12 7T 3 



d 2 q 



K^°(k,q) 



:i-e 



(fe — q) 2 q 2 



,2 A 



(73) 



where we have abbreviated the action of the NLO kernel on the left hand side of Eq. (1721) . To 
perform the part of g-integral involving Kj LO (k,q) we will employ Eq. (164)) . The rest of the 
g-integral is 



2 k 2 

d 2 q Tl r~ In — q 

[k — q) z q z 



2 A 



2vr [ip'(l + A) - V'(-A)] jfe 



2 A 



(74) 



which can be derived using Eq. (1191) and the formulae in Appendix [A] Combining this with 
Eq. (JMD yields 



a n K f 



X 



-2 A 
01 



o _ 2A+1 r(i-A) rfk 
r(A) y fc 2 (i 



3 'ifc-£C01 > 



a 2 N Nf f 1 

"M iVciV / fc 2A J 1 
12 7T 3 



X 2 (-A) + x(-A) In ^ - i ^(-A) + J + A) 



1 

2 



(75) 



Performing the integration over A; we obtain 



x 



-2A 



01 



\- 2 x(-x)+x(-x) 



4 e ~ 2 7 2 

ln — — 2 -x(-A) - - 

. ^01^ A 



-^'(-A) + ^'(l + A)S>x - 1 2A , 



(76) 
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where 7 = — is the Euler's constant. Since the LO BFKL kernel acting on a power in 
transverse coordinate space gives the usual BFKL eigenvalue 



we write 



x 



-2A 
01 



7i 



X(-A)x 



-2 A 
01 



(77) 



LO 



^NLO] 



X 



-2 A 
01 



7T 



X(-A) 



a, 



.ft In 



4e 



-27 



X 0i/4ts 



127T 



A B (A) 



-2 A 
01 



(78 



with 



A B (A) 



V 



-A) + ^ X (-A)+^(-A)-^(l + 



A) + yX(-A). 



(79) 



Here one might worry that this value of A#(A) depends on our choice of the constant under 
the logarithm in Eq. flTHj) . However this choice is not arbitrary and is consistent with the 
constant obtained under the logarithm of the coordinate space running coupling corrections 
resummed to all orders in [27]. (In [27] we would also have a factor of e _5//3 under the running 
coupling logarithm. Here, following the convention of [29,30,34] we have chosen to place this 
contribution separately: it leads to the y x( — A) term in Eq. (1791 .) 

Replacing A = 7 — 1 reduces Eq. (1791 to Eq. (1691 ! We have thus shown that our intercept 
(1680 and the intercept in Eq. (1691 found in [29] are consistent with each other. The difference 
between them is due to the fact that they are calculated for two different observables. 



8 Conclusions 

In this paper we have studied the implications for the linear BFKL evolution equation of 
including running coupling corrections into the non-linear JIMWLK and BK evolution equations 
performed in [27]. In particular we have derived the BFKL equation including running coupling 
corrections to all orders and calculated the leading-iV/- NLO BFKL pomeron intercept. Until 
now the form of the BFKL equation with all orders resummed running coupling corrections 
was not known. What existed was a conjecture by Braun and by Levin [39,40] based on an 
interesting assumption that bootstrap equations hold even with running coupling corrections 
included. We have shown that this conjecture is in fact accurate, though for a slightly different 
observable than suggested originally. Our results were derived using the s-channel language of 
light cone perturbation theory [47, 48] . Up to now the NLO BFKL intercept had only been 
calculated either by using the standard Feynman perturbation theory in [30, 34] and or by 
employing the background field method [29] . Our calculation provides an independent check of 
the intercept found in [30,34] and connects it to the one obtained in [29]. 

Our paper has two main results. First of all, we have obtained the BFKL equation for the 
unintegrated gluon distribution including running coupling corrections resummed to all orders. 
The equation is given by (1461) . After a redefinition of the unintegrated gluon distribution shown 
in Eq. (j4"7j) . Eq. (j4*6|) leads to Eq. (j4*8l . which was conjectured by Braun and by Levin in [39,40] 
by postulating the bootstrap condition in the running coupling case. We have thus shown that 
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the conjecture of [39,40] for the BFKL equation with the running coupling corrections is correct, 
though for a slightly non-traditional definition of the unintegrated gluon distribution (14 7p . We 
have clarified that this result is based on a complete absorption of all a s Nf corrections into the 
running of the coupling in the BFKL equation (as in Eq. (|46p ). This is one of many possible 
constructions as explained in Sec. [31 all others would separate off explicit A^- contributions that 
are not interpreted as running coupling contributions. The result Eq. (I46p is formulated in the 
form of a "triumvirate" of the couplings [27]. 

Our second result is an independent check of the leading-A/ contribution to the NLO 
BFKL intercept (expanded to strictly NLO, without any resummations), which we performed 
in Sect. El Our intercept, obtained for the evolution of the unintegrated gluon distribution 
function, is given by Eq. (|68|) and completely agrees with the results of Camici and Ciafaloni [34] 
and of Fadin and Lipatov [30] . The NLO intercept appears to strongly depend on the physical 
observable: we demonstrate that by calculating the leading-A/ NLO BFKL intercept for the 
dipole scattering amplitude, starting from that of the unintegrated gluon distribution. Our 
result is given in Eq. f!79l) and agrees with the result of Balitsky [29]. The difference in the 
expressions fl67j) and f!79|) is fully explained by the fact that the two intercepts refer to different 
observables. 
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A Useful formulae 

Here we list some useful mathematical formulae. The main formula we need above is 



oc 




(Al) 







Using Eq. flAll) one can derive the following useful results (here b > 0): 
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B Evaluating Eq. (EH) 



Our goal here is to perform the fc-integration in Eq. (l27j) . We begin with the first term in the 
curly brackets in Eq. (l27j) . Performing the a-integral first we obtain 



d\_H_fk_ o _ iq . (zi _ aio)+iq . (zi _ aJl) 1 (q-k).(qi-k ) (q - fc) 2 



7 ' (2tt) 2 (2tt) 2 (27r) d 
Eq. (1B1I) can be rewritten as 
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Defining a new integration variable k = k — (1 — (3) q — /3 q' we get 

i 



t2 j2 / 

A N* I q q p -ilizx-x a )+iq-{zx-xx) 

f 1 (27r) 2 (27r) 2 



q 2 q' 2 



d d k k 2 -f3(i-f3)(q- q r (m 

(27r) d ~k 2 + P(l-P)(q-q<) 2 ' 



where we dropped the terms linear in k in the numerator as they vanish after angular integra- 
tion. Performing the fc-integral yields 
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d2( l rf V p-iq-izx-^+iqizx-nx) 



(2tt) 2 (2tt) 2 



4* / 



d/2 



(B4) 



Inserting d — 2 — e, expanding the expression in powers of e, replacing 1/e with ln/^fs and 
integrating over /? yields 
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where fi 2 -^ = fi 2 MS Air e~ 7 . Eq. (1B5[) gives us the first term in the curly brackets of Eq. (|28|) . 

The last term in the curly brackets of Eq. (1271) gives us zero after performing a dimensionally 
regularized fc-integral. We are left only with the second and the third terms in the curly brackets 
of Eq. (12?]) . which are analogous to each other. Here we will show how to do the fc-integration 
in the second term only: the integral in the third term can be easily done in the the same way. 

Defining k = k — (1 — a) q the second term in the curly brackets of Eq. (1271) can be written 

as 

- 4N , fi«iv t *-*,-,> r da ( B6) 



(27r) 2 (2vr) 2 q 2 q' 2 J J{2^) d k^ + a(l 
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which, repeating the above steps which led from Eq. (1B3[) to Eq. (1B5p . can be recast into 

N f f d 2 q d 2 q> e _ iq . {zi _ Xo)+iq<Zl _ Xl) J_^fe^^ ^ 



3vr J (27r) 2 (27r) 2 q' 2 fi 2 



MS 



This is exactly the second term in the curly brackets of Eq. (128]) . The third term is done by 
analogy with the second one. 
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